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We derive an index theorem for zero-energy Majorana fermion modes in a superconductor- 
^v^j \ topological insulator system in both two and three dimensions, which is valid for models with 

chiral symmetry as well as particle-hole symmetry. For more generic models without chiral 
symmetry, we suggest that Majorana zero modes are classified by Z2. 
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Majorana bound states have recently been attracting much interest as elementary ex- 
citations obeying non-Abelian statistics.® They have been predicted in v = 5/2 quantum 
Hall effects,® p-wave super conductor sJH3 topological insulators with the proximity effect of 
s-wave superconductors® as well as of ferromagnetic insulators, 7 ' s-wave superfluids of ul- 
tracold fermionic atoms,® vacuum states of superfluid 3 He-B,® superconducting states in 
graphene,®^ and massive Dirac equations in magnetic fieldsP® Recently, Teo and Kane^ 1 
have proposed a model of superconductor-topological insulator systems in three dimensions. 
They have shown that the third Chern number, which is a topological invariant for three- 
dimensional topological insulators,®^ 1 reduces to a winding number of order parameters on 
a two-dimensional sphere surrounding a pointlike topological defect. This suggests that the 
. Majorana zero modes have a topological origin and are stably protected from perturbations. 

This property may play a crucial role in the application of Majorana zero modes to topological 
quantum computations p3ED> Along this line, Tewari et. al. have established a theorem for the 
existence of Majorana zero modes in chiral p-wave superconductors 51 and in spin-orbit-coupled 
semiconductors 

^ ■ In this paper, we derive an index theorem for the Majorana zero modes in a model of 

superconductor-topological insulator systems. To be concrete, we establish the relationship 
between the sum of chiralities of the Majorana zero modes and a topological invariant asso- 
ciated with order parameters. This indicates that the Majorana zero modes in such systems 
are stable owing to their topological origin. When we derive the above index theorem, we use 
a simplified model that has not only particle-hole symmetry, as the Bogoliubov-de Gennes 
Hamiltonians should have, but also an additional chiral symmetry. We next discuss the effects 
of chiral symmetry breaking perturbations, especially those of the chemical potential term. 
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We suggest that the number of zero modes is classified by Z2. Namely, only the evenness or 
oddness of the number of zero modes in models with chiral symmetry is invariant under chiral 
symmetry breaking perturbations. This is consistent with the results obtained by Tewari et. 
o/P Finally, we present a concrete zero mode wave function for the Fu-Kane modeP* in both 
the cases of fi = and fi 7^ 0. 

Let % be a Hamiltonian defined by 

n = w = i^d j + 1 a ( f> a , (i) 

where j = 1, • • • , d and a = d + 1, • • • , 2d, and 7 M = 1, • ■ ■ , 2d) are 7 matrices in 2d di- 
mensions {7^,7^} = 25^ u . 4> stands for a set of order parameters. This Hamiltonian has been 
introduced by Kane and his collaborators^^ as a minimal model describing superconductor- 
topological insulator systems. In the case of d = 2, cj> = Ao(r)(cos 6>, sin0) is the model 
proposed by Fu and Kane 6 ' to study a Majorana zero mode in a vortex of an s-wave super- 
conductor, where (r, 9) are the polar coordinates in two dimensions. In the case of d = 3, 
4> = (ReA,ImA,m) describes generic superconductor-topological insulator systems studied 
by Teo and Kanep^ 1 where the region m < corresponds to a bulk topological insulator, with 
the proximity effect induced by a superconductor. If 7^ is constant, the spectrum of 

the Hamiltonian has a gap, describing a massive Dirac particle. Even with a more generic <f> 
that is not spatially uniform, the spectrum may still be massive. However, if 4> has a topolog- 
ical defect, there appear bound states around it. In particular, it has been showrP^ that the 
present model (pQ) can admit generically bound states sitting exactly at zero energy. To study 
such zero modes bound to a pointlike defect, we assume that away from the defect located at 
x = 0, the asymptotic form of the order parameter (f) is given by 

101 4>o = const. (|x| — > 00). (2) 

The explicit expression for 7 matrices is not necessary,^ but we have to arrange them to 
obey Cj^C~ 1 = 7^ and C~f a C~ l = —7°, where C (C 2 = 1) is the anii-unitary operator^ 
describing particle-hole symmetry. Therefore, the Hamiltonian ([T]) has particle-hole symmetry, 

CUC- 1 = -H, (3) 

as it should be in the case of the Bogoliubov-de Gennes Hamiltonians. The wave function 
of the zero modes can be chosen such that ipQ t i(x) = Ctpoj(x), which has been referred to 
as a Majorana fermion mode. Another symmetry of the Hamiltonian ([1]) is chiral symmetry 
denoted by 

75^75 = -n, (4) 

where 75 = (—i) d ^f 1 ■ ■ -j 2d with the property C^sC -1 = 75. It should be noted that the 75 
matrix is that in 2d dimensions in the ordinary sense, although we are dealing with a model 
in d spatial dimensions. This symmetry seems artificial, since the chemical potential term and 
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many other terms with some products of 7 matrices break chiral symmetry. According to the 
Altland-Zirnbauer classification schemep^^ the model only with particle-hole symmetry (|3|) 
is classified as class D, which is the most generic class of superconductors, whereas the model 
with an additional symmetry @ belongs to class BDI, which is the time reversal invariant 
chiral class PED 

Even with such an enhanced symmetry (and although we have neglected background gauge 
fields for simplicity) , the minimal model ([I]) is nevertheless useful for clarifying the topological 
stability of the Majorana zero modes, as we shall show momentarily. When d = 2, it has a 
deep relationship with a Dirac fermion coupled to a Higgs field with a vortex background 
Jackiw and Rossi 22 ' have derived the zero mode wave functions of this model explicitly. It has 
been shown by Weinberg®) that these zero modes are indeed ensured by the index theorem, 
even though this model is defined for the noncompact two-dimensional space R 2 . The model 
in three dimensions is associated with a Dirac fermion coupled to a Higgs field with a 't 
Hooft-Polyakov monopole^"^ background whose index theorem, surprisingly valid for open 
and odd-dimensional spaces, has been derived by Callias. 26 ' In what follows, we first prove 
the index theorem for the Hamiltonian ([1]) simultaneously in d = 2 and d = 3 cases using the 
method developed by Weinberg P^ 1 and next, we discuss the effect of the chemical potential 
that breaks chiral symmetry. 

Let us first define the index of the Hamiltonian (pQ) by ind% = N+ — N-, where N± denote, 
respectively, the number of zero modes of ([1]), = 0, with ± chirality, Ts^o^ = ±^oi- 

This can be written as 

2 

Tfl 

ind'H = lim Tr75 (5) 

m->-0 —u A + m A 

where Tr denotes the trace over the space coordinates as well as the 7 matrices. On the other 
hand, it is well-known that in the large mass limit, the same quantity yields the topological 
invariant 

M 2 

Crf = iToo Tr75 -^ + M2 ' (6) 
where Cd is the dth. Chern number^ 1 if d = even, while Cd = if d = odd. It is possible to 
express Tr in both eqs. ([5]) and ([6]) using the eigenstates of the Hamiltonian, l~Lip n — X n ip n . 
In the case of nonzero energies, <^_„ = 75^ has the energy A_ n = — X n , and hence, these 
two states are orthogonal to each other, J d d x(fn^ip n = J d d xf\ l ip- n = 0. Then, it is easy to 
see that eqs. ([5]) and ([6]) give the same index of the Hamiltonian. However, in the plane-wave 
basis, the above two quantities become two kinds of different invariants, giving a nontrivial 
relationship between them. The necessity of two kinds of masses is manifested when the axial- 
vector current is introduced such that 

J i (x,m,M) = lim tr7 5 7 i (— — — ) S(x - y) 

y-+x \L> + m D + M J 
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where tr denotes the trace over the 7 matrices, and in deriving the second line, we have used 
the fact that T> anticommutes with 75. The first term is, in the limit of m — > 0, the axial- vector 
current of the model under consideration, while the second term is the Pauli-Villars regulator 
whose mass M is set to infinity after the calculation. Such a regulator field is needed to make 
the current well-defined. 

It is well-known that the index ([5]) and the Chern number §§§ have an intimate relationship 
with the chiral anomaly, which yields an anomalous term to the conservation law of the 
axial-vector current. To see this, it is convenient to describe the current on the basis of 
eigenstates of the Hamiltonian. The completeness of eigenstates of the hermitian Hamiltonian, 
Y.n <Pn(x)<Pn(y) = S(x - y), leads to 

J l {x, m,M) = -J2 ^ n {xh 5 f(-i\ n ) ( A 2 + m 2 ~ A 2 +m 2 ) {Pn ^' 

Then, the divergence of the current is easy to compute, since the derivative di operates only 
on the wave functions in the above expression. After some calculations, in particular, using 
the relation j J dj(p n = i{^ a 4>a — A n )</3 n , we obtain 



It thus turns out that the divergence of the current (J7J) is given by 

( m 2 M 2 \ 

d^(x,m,M) = lim(-2)tr 75 [-^^ ~ _ p2 + M , ) ^ ~ v)- 

Integrating both sides over the space R d and taking the limits m — > and M — > 00 yield 

mdU-c d = -^J d d xd i J i {x^,oQ) 



dSiJ\x, 0,oo), (8) 



where dS{ is the infinitesimal surface element on S rf_1 , i.e., the boundary of R d at |x| — > 00. 
This equation indicates that in even dimensions if the space manifold is closed, the surface 
term above vanishes, and the index of T~L is just given by the Chern number c^. For generic 
open manifolds, the index theorem is modified by the surface term given directly by the axial- 
vector current. It should be noted that in the present model ([T|), Cd always vanishes since we 
have neglected the gauge fields, for simplicity, and thus, the index of the Hamiltonian, which 
is due to the Majorana zero modes, is given not by the Chern number, but by the surface term 
in eq. ([8|). This surface term is also a topological number associated with the configuration of 
the order parameter (f> at \x\ — > 00, as we shall see below. 

To compute the surface term, we now calculate the axial-vector current ([7]) in the plane- 
wave basis. To this end, note — V 2 = —d 2 + 4> 2 a + j^ a (idj4> a )- Therefore, at \x\ — > 00, we 
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have 

e - ikx (-V 2 )e ikx = k) + eg + 7-V(^a) - d] - 2ik j dj, 
where we have used eq. ([2]). This leads to 
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tr 75 y e - ikx T> 



d d k 
(2vr) d 



1 



1 



-V 2 + m 2 -£> 2 + M 2 



ikx 



k 2 + g + m 2 + -ih a {idj(t>a) + k 2 + 4> 2 + M 2 + j^ a {id j( /) a ) + O 



/^75f(^ + ^)E( 



(fc 2 + g + m 2 )^ 1 (A; 2 + 2 + M 2 )«+! 



x [ 7 V(^ a ) + 0] 



where O = —d 2 — 2ik 3 dj. Note that tr 7 5i, where T is a product of 7 matrices T = 7 ^ 7 y • • • , 
can be nonzero only when T includes at least as many as 2d 7 matrices, which occurs in the 
n = 1 (n = 2) term for d = 2 {d = 3) in the above expansion. Using tr 7 5 7 /il 7 ^ 1 • • • r ) lLd ^ Vd = 



{2ife 



, we have in the case of d = 2, 

■7^,0,00) = * tr 757 V*7>a) 7 V(^A) + 0(M" 2 ) 



-e^e ab ^> fe + 0(|x|- 2 ) 

7T 



where 



t> a /(f)Q is the normalized field, 



1, at |x| — > oo. Note here that dj</> a ~ 0(|x| x ^ 



since |0| approaches a constant value. In the case of d = 3, 

^-3tr 757 i (-i 7 > a )y 7 6 (ia i 6 ) 7 fe 7 c (^ fc c ) + o(|x|- 3 ) 



J J (x, 0, oo) 



32vr< 



'o 



j-e^ k e abc 4> a d^ b dj c + C>( 



Substituting these into eq. 0, we obtain 

( 1 



ind^ 



^ / dx'^^di 
1 

8^ 



dSie ijk e abc 4> a d^ b d k ^ c 



(d = 2) 
(d = 3) 



(9) 



Since in eq. Q, <^ is a unit vector with d components defined at \x\ — > oo of R d , (j> can be 
regarded as a mapping from S ' -1 to S d_1 . It is well-known in mathematics that iTd-i(S d ~ 1 )=Z, 
implying that such a mapping is classified by integers. The r.h.s. of eq. Q denotes exactly 
the winding number of </> over S^ 1 at |x| — > oo. 

In the case of the Pu-Kane model,® we can choose = (cos q9, sin qO) at r — > oo if we 
consider a generic vortex with vorticity q. Substituting this into eq. Q, we then see that 
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the winding number in eq. © is exactly —q, implying that the Hamiltonian has at least 
q zero modes. It is not difficult to solve the eigenvalue equation and to show that such a 
vortex indeed admits exactly q zero modes. The r.h.s. of eq. in the case of d = 3 has been 
derived by Teo and Kane^ as a topological invariant for the many-body ground state of the 
model. They have calculated the third Chern number^ associated with non-Abelian Berry's 
connection and shown that it reduces to the winding number of the order parameter over S 2 . In 
the present calculations, without resorting to adiabatic assumptions, we have established the 
index theorem that is directly related to the analytic invariant associated with the Major ana 
zero modes with the topological invariant of the order parameter. 

Thus far, we have derived the index theorem for a superconductor-topological insulator 
model in the case where it has not only particle-hole symmetry but also chiral symmetry. 
However, let us consider the term Ti' = i"f a "f b n a b, for example. We see that it transforms as 
C7~L'C~ l = —H', but 75%'75 = +%' ■ Therefore, generic Bogoliubov-de Gennes Hamiltonians 
have broken chiral symmetry. In particular, the chemical potential is one of the most important 
terms breaking chiral symmetry when we study the superconductivity. Hence, we will next 
investigate the effect of the chemical potential term — /i/3 that breaks chiral symmetry, where 
(3 is a matrix with the transformation properties = — (3 and 7s/37^" 1 = +/?P^ When 

\x is a small constant, the perturbation theory indicates that the leading order corrections of 
the energies are given by the eigenvalues of the matrix H 1 ^ — — fi{ifo^\(3\ipQj) , where |y?o,i) 
denotes the ket notation of the unperturbed ith zero mode wave function. In general, we can 
choose l^oj) = C\ipoj), as already discussed. Then, we see that each matrix element is purely 
imaginary, %[*■ = —H^j because of the particle-hole symmetry of (3. It thus turns out that H' is 
an antisymmetric purely imaginary matrix. Without further symmetries, off-diagonal matrix 
elements of Ti' do no vanish genericallyPS This leads to an interesting suggestion® that if an 
unperturbed (/x = 0) system admits several zero modes, some of them couple together through 
the chemical potential term and become nonzero energy states such that two states form a 
pair with opposite energies. This is due to particle-hole symmetry. Namely, an even number of 
zero modes become nonzero energy modes if the chemical potential term is switched on. This 
implies that in a system with an odd number of unperturbed zero modes, there remains at 
least one zero mode even if any perturbations with particle-hole symmetry are added. In other 
words, the rank of an antisymmetric matrix is always even, implying that odd-dimensional 
%' has inevitably an odd number of zero eigenvalues. We thus conclude that Majorana zero 
modes bound to a pointlike defect are classified by Z for class BDI, whereas they are classified 
by Z2 for class D. This is consistent with the results obtained by Tewari et. a/. 5 ' 

Finally, we present, as a nontrivial example, an explicit wave function of a zero mode in 
the Fu-Kane modeP 1 with a q = 1 vortex. With the 7 matrices given inp» the zero mode 
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wave function obeys 

/ iaidj-fi <f>i-i(f>2 \f S)_ Q 

\ <pi + i<p 2 -ia j dj + n J \ i] J 

Note that <pi — icf>2 = Ao(r)e~ tS , and that particle-hole symmetry makes it possible to choose 
r\ = —<T2^*^ Then, the above equation reduces to 

(ia j dj A (r)e- ie a 2 e = 0. (10) 

If we further set £ T = (v, u), we see that when /j, = 0, two equations for u and v are decoupled, 
and the normalizable solution is given by v = and 

u(r,8) oce-f rdr ' A °( r '\ 

which was obtained in RefP Here, we have assumed that Ao(r) > for r — > 00. This wave 
function has chirality —1, and hence, the index of this model is —1. The above discussions on 
the classification scheme for classes BDI and D suggest that the present zero mode is stable 
against perturbations breaking chiral symmetry. Indeed, even when /i 7^ 0, eq. (|10p has the 
following normalizable wave function: 

\ Jo(j*r) J 

where J n (x) is the Bessel function. 

In conclusion, we have derived the index theorem ([9]) for the model of a superconductor- 
topological insulator system in class BDI proposed by Kane et al. This ensures the topological 
stability of the Majorana zero modes bound to a topological defect. For a model in class D, i.e., 
a model with particle-hole symmetry only, we have suggested that zero modes are classified 
by Zg. Namely, only the parity (evenness or oddness) of the number of zero modes (the index, 
in other words) in class BDI is invariant against chiral symmetry breaking perturbations. 

We would like to thank A. Furusaki for useful comments on the universality classes of 
topological insulators and superconductors, and also C. L. Kane and the Referee of JPSJ for 
pointing out the Z2 classification of zero modes in class D. This work was supported in part 
by Grants-in-Aid for Scientific Research (Nos. 20340098 and 21540378). 



7/1 



J. Phys. Soc. Jpn. Letter 
References 



1 
2 
3 
4 
■5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 



N. Read and D. Green: Phys. Rev. B 61 (2000) 10267. 
D. A. Ivanov: Phys. Rev. Lett. 86 (2001) 268. 

A. Stern, F. von Oppen, and E. Mariani: Phys. Rev. B 70 (2004) 205338. 

S. Das Sarma, C. Nayak, and S. Tewari: Phys. Rev. B 73 (2006) 220502. 

S. Tewari, S. Das Sarma, and D.-H. Lee: Phys. Rev. Lett. 99 (2007) 037001. 

L. Fu and C. L. Kane: Phys. Rev. Lett. 100 (2008) 096407. 

Y. Tanaka, T. Yokoyama, and N. Nagaosa: Phys. Rev. Lett. 103 (2009) 107002. 

M. Sato, Y. Takahashi, and S. Fujimoto: Phys. Rev. Lett. 103 (2009) 020401. 

G. E. Volovik: Pis'ma ZhETF 90 (2009) 639. (iarXiv:0909.3084|l 

I. F. Herbut: arXiv:0909.4231. 



I. F. Herbut: larXiv:0910.4906l 
J. C. Y. Teo and C. L. Kane: arXiv:09094741. 

X.-L. Qi, T. L. Hughes, and S.-C. Zhang: Phys. Rev. B 78 (2008) 195424. 
A. Kitaev: Ann. Phys. 303 (2003) 2. 

C. Nayak, S. H. Simon, A. Stern, M. Freedman, and S. Das Sarma: Rev. Mod. Phys. 80 (2008) 
1083. 

16) S. Tewari, J. D. Sau, and A. Das Sarma: larXiv:0910.4763l 

17) In the Fu-Kane model, they have used 7 1 = a 1 ® a 3 , -f 2 = a 2 ® a 3 , -y 3 — 1 (g) cr 1 , and j 3 = 1 <g) a 2 , 
where in the notation of a ® 6, a and b denote the spin and particle-hole spaces, respectively. 
The particle-hole transformation C is defined as C = a 2 ® a 2 K, where K denotes the complex 
conjugation operator. In the Teo-Kane model, 7 J = a 3 ® a 3 ® a 3 (j — 1, 2, 3), 7 4 = 1 <E> a 1 ® 1, 
7 5 = 1 ® a 2 ® 1, and -f 6 = 1 (g> a 3 ® cr 1 , where in a ® b ® c, a, 6, and c denote, respectively, the 
spin, particle-hole, and orbital spaces. For this model, C — a 2 (g> cr 2 ® IX. 

18) M. Zirnbauer: J. Math. Phys. 37 (1996) 4986. 

19) A. Altland and M. Zirnbauer: Phys. Rev. B 55 (1997) 1142. 

20) A. P. Schnyder, S. Ryu, A. Furusaki, and A. W. W. Ludwig: Phys. Rev. B 78 (2008) 195125; 
larXiv:0905.2029l 

21) A. Kitaev: larX"iv:0901. 26861 

22) R. Jackiw and P. Rossi: Nucl. Phys. 190 (1981) 681. 

23) E. J. Weinberg: Phys. Rev. D 24 (1981) 2669. 

24) G. 't Hooft: Nucl. Phys. B 79 (1974) 276. 

25) A. M. Polyakov: JETP Lett. 20 (1974) 194. 

26) C. Callias: Commun. Math. Phys. 62 (1978) 213. 

27) For a review, see, e.g., K. Fujikawa and H. Suzuki: Path Integrals and Quantum Anomalies (Oxford 
University Press, Oxford, 2004). 

28) In the Fu-Kane model, /3 = -iy 3 ^. 

29) If the matrix /3 has chiral symmetry, H 1 vanishes in general. This alternatively indicates the 
stability of the index ©. 

30) This condition is represented as tp c = Ctp = —if, where (p T — {£, T ,r) T ). The factor —i is for 
notational convenience. 



8/0 



